Quantum Gaussian Channels with Additive Correlated Classical Noise 
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We provide a model to study memory effects in quantum Gaussian channels with additive classical 
noise over an arbitrary number of uses. The correlation among different uses is introduced by 
contiguous two-mode interactions. Numerical results for few modes are presented. They confirm 
the possibility to enhance the classical information rate with the aid of entangled inputs, and show 
a likely asymptotic behavior that should lead to the full capacity of the channel. 
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I. INTRODUCTION 

A memoryless quantum communication channel makes 
the fundamental assumption that the noise between con- 
secutive uses of the channel is independent. In many 
real- world applications this assumption may be good, but 
for many others the noise may be strongly correlated be- 
tween uses of the channel. Such a possibility has been 
recently put forward in channels with continuous alpha- 
bet, specifically lossy bosonic channels 0. The main 
motivation that has led to investigate memory effects in 
such channels has been the possibility to enhance their 
classical capacity by means of entangled inputs 0. 

Lossy bosonic channels belong to the class of Gaussian 
channels which play a prominent role because it might be 
simple enough to be analytically tractable, thus providing 
a fertile testing ground for the general theory of quantum 
capacities |3|. Another example of Gaussian channel is 
provided by a channel that adds thermal noise to the 
input. For that channel it has been recently proved that 
entangled inputs enhance the 2-shot classical capacity in 
presence of correlated noise |4( . 

Here we present an extension of the model used in 
Ref. , that can be employed for an arbitrary number of 
channel uses. The memory effect is modeled by assum- 
ing that the noise affecting subsequent uses of the chan- 
nel follows a Gaussian distribution with correlation in- 
troduced by contiguous two-mode interactions. We show 
that if the memory is non-zero and the input energy con- 
straint is satisfied, the channel transmission rate can be 
improved by entangled input states instead of product 
states. Moreover, we notice the appearance of an asymp- 
totic value of the transmission rate after many uses of 
the channel that should lead to the full capacity of the 
channel. 

The paper is organized as follows. In Sect. II and 
III, we review the basic notions about bosonic channels 
and Gaussian channels with additive classical noise. In 
Sect. IV we describe the action of the memory. Finally, in 



•Electronic address: ruggeri@le.infn.it 
^Electronic address: stefano. manciniOi 



Sect.V we analyze the transmission rate of this channel 
and discuss our results. 



II. BOSONIC CHANNELS 

A bosonic channel is any completely positive and trace 
preserving map acting on the state of a mode of electro- 
magnetic field. On multiple uses of that channel one has 
to account for many of such modes. Let us consider n of 
them associated with n pairs of annihilation and creation 
operators aj and Oj (satisfying the canonical commuta- 
tion relation), or equivalently of quadrature components 



q.j = (aj + a])/\/2 and pj 
these in a column vector 



a\)/y/2. Ordering 



r = [qi 



,q n ,pi, 



,Pn\ 



(1) 



we may define the mean vector m and the covariancc 
matrix V of an n-mode state p as 



V = Tr 



m = Tr (pr) , 



(r — m) p (r — m) 



(2) 



(3) 



where J = i f ^ ^ is the symplectic matrix, with X 

the n x n identity matrix. 

If the state p is Gaussian is completely characterized by 
the mean vector m and the covariance matrix V. More- 
over, its von Neumann entropy is given by 



s( P ) = J2g(\x j \-l), 



(4) 



where ±Aj are the symplectic eigenvalues of the covari- 
ance matrix V, that is the solutions of the equation 



and 



unicam.it 



det [V - \J] = 0, (5) 
(x + 1) log 2 (x + 1) — x log 2 x x > 







x = 



(6) 



2 



is the entropy of a thermal state with a mean photon 
number x. 

Let us now consider the classical use of a bosonic chan- 
nel described by a completely positive and trace preserv- 
ing map T: 



P >-> T[p], 



(7) 



Then, quantum states carry the values of a random clas- 
sical variable. Quite generally we could think to map 
phase space points into quantum states, hence for each 
mode we consider two real random values (or one com- 
plex). We then label an input state over n modes (uses) 
as pa with a. E C™. 

The coding theorem for quantum channels asserts that 
the classical capacity of the bosonic channel T is given 
by 



C (T) = lira — max 



S^J P(d 2n a)T[p a 
J P(d 2n a )S(T[p a }) 



(8) 



where S denotes the von Neumann entropy and P(d 2n a) 
is a probability measure. 

In Eq.©, the maximum is taken over all probability 
measures {P{d 2n a)} and collect ions of density operators 
{pa} satisfying the energy constraint 

i J P(d 2n a)Tr (p a J2a]aA <n, (9) 

with n the maximum mean photon number per mode at 
the input of the channel. 



III. GAUSSIAN CHANNELS WITH ADDITIVE 
CLASSICAL NOISE 

Consider now a specific bosonic channel T acting as 
follows 



T[p a ] = J d 2n f3 Q((3) D((3)p%DH(3), 



(10) 



where D((3) = £»(/3i)£»(/3 2 ) ■ ■ • D((3 n ) with D[fij) denot- 
ing the displacement operator on the jth mode, (3j G C. 
This is known as Gaussian channel with additive classical 
noise once the kernel is Gaussian, e.g. 



q{p) = n 



fe=i 



i 



(11) 



In such a case the channel randomly displaces each input 
state according to a Gaussian distribution, which results 
in a thermal state (N is the variance of the added noise 
on the quadrature components qj and Pj, or, equivalently 
the number of thermal photons added per mode by the 
channel) . 



Since for Gaussian channels it is conjectured that 
Gaussian inputs allow to achieve the capacity Q , let us 
consider to encode a into coherent input states 



p% = D{ a )\0)(0\rt( a ), 



(12) 



where D(a) = D(ai)D(a 2 ) . . . D{a n ) and |0) 
f) . . . |0) are the vacuum states of the n modes. More- 
over, suppose that the state (|12f) is drawn with probabil- 
ity 



n 1 

p(a) = n — 



k=l 

then the output states and their average read 
f& = T |>S] , 
-p° ut = J d 2n oc P{ a )p™\ 



(13) 

(14) 
(15) 



The Gaussian map effected by the Gaussian channel 
(|10fl . in case of Gaussian input state, can be solely char- 
acterized by the covariance matrices, that is 



V 



V m + diag(AT,iV,...,A0, 



(16) 



where V m is the input covariance matrix. For the states 
of Ea.iPJl. it is 



V m = idiag(l,l,. 



(17) 



We could also consider the possibility of entangled in- 
put states. These can be accounted for by the transfor- 
mation pJ5 — > S(7?.)pqE' 1 ' (1Z) that uses the multimode 
squeezing operator [5| 



E (R) = cxp 



1 - 



. > ^ y n ika]al- K* jk >i,,n. 
jk=i,j^k 



, (18) 



with 1Z a symmetric n x x matrix. Thus, the input co- 
variance matrix results 



V 



1 / cxp (n) o 

2 1 exp(-Te) 



(19) 



Let n r be the average number per mode of squeezed pho- 
tons added at input, i.e. 



n r 



n 

-$>|E(ft) a J aj Et(ft)|0). 



(20) 



3=1 



Then, due to the bound the effective number of pho- 
tons to input on each mode would be (n — fir). 



IV. THE MEMORY MODEL 

Now, let us consider memory effects in the Gaussian 
channel with additive classical noise. The output state 



3 



(|10f> can be seen as the convolution of the input Gaussian 
state l|12[l with a Gaussian thermal state having average 
number of thermal photons N per mode. To introduce 
correlations in the latter we use arguments from mul- 
timode squeezing. It is physically reasonable that only 
contiguous uses of the channel would be strongly corre- 
lated, thus we introduce the following symmetric matrix 



S 



( 1 

1 1 

■•101 ■•■ 
1 1 

V o i o / 



(21) 



where s is the parameter measuring the degree of mem- 
ory. Then, we take the added noise as characterized by 
the covariance matrix 



V 



vi 



N 



where V-f is the diagonal matrix 



[v,\ = Jv- £ [v 2 \, 



and 



V * = i 



exp (S) 
exp (-S) 



(22) 



(23) 



(24) 



In order to realize a physical transformation, the channel 
must have V N > 0. This is clearly satisfied if > 
and > 0. The first condition determines the range 
of allowed values for the memory parameter s, while the 
second condition is always satisfied. The presence of the 
parameter e in Eq. l|22|l only deserves as mathematical 
trick to guarantee the positivity of V N even in the case 
of N < 1/2. It is 



e = 1 
< e < 2N 



N>1 
~ 2 



(25) 
(26) 



Notice that the noise correlations so introduced are clas- 
sical and for s — ► in Eq. (|22|l we recover the memoryless 
case. 

For the entangled inputs we consider in Ea. (|18fl the 
symmetric matrix 



K = 



( ° 1 
1 

... i 



;■ : ■: \ 

1 ••• 

1 1 

V o i o ) 



(27) 



with r is the input entanglement parameter. 

At the output of the channel, we get states with the 
covariance matrix 



R 
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FIG. 1: Transmission rate R (maximized over the correlation 
coefficient y) as a function of the input entanglement param- 
eter r for a channel with a degree of memory s = (dotted 
lines), s — 0.1 (dashed lines), s — 0.2 (solid lines). The num- 
ber of uses is from top to bottom n — 2,3,4,5. The input 
mean number of photons is n = 2, while the number of added 
photons is N = I. 



In turns, the covariance matrix associated with the mix- 
ture of the output states lfl"5|) would be 



V" 



V 1 



V 



N 



(28) 



v = v° 



(29) 
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Here 

K = K 1 +9K i . (30) 

with JC\ a diagonal matrix 

[Ki]jj = (n - n r ) - 9[K 2 ]jj, (31) 

and K.2 being of the same form of with s — ► —y. The 
coefficient y accounts for possible classical input correla- 
tions. 

Also in this case we must have JC > 0. This is clearly 
satisfied if JC% > and IC2 > 0. The first condition de- 
termines the range of allowed values for the parameter 
y, while the second condition is always satisfied. Like- 
wise the parameter e in Ea. (|22[l . we have introduced in 
Ea. (|30f) a parameter 9 that only deserves to guarantee 
the positivity of JC when (n — n r ) < 1/2. ft is 



9 = 1 

< 9 < 2(n-n r ) 



-s 1 

(n-n r ) > - 

r- - \ 1 



(32) 
(33) 



According to Eq.JSJ), the capacity is given by the 
asymptotical behavior (n — + 00) of the maximum of 
transmission rate 



R 



1 



[S (p out ) - S (p out )] 



(34) 



Since p out and p out are both Gaussian states, following 
Eq.Q we have 



R(r,y) = - 



(35) 

Thus, in order to evaluate the transmission rate, we need 
to compute the n symplectic values Xj , \° ut of V , 
V out respectively. 



RESULTS AND CONCLUSIONS 



The transmission rate of Ea. (|35|l has been evaluated 
numerically restricting to some cases where e — 9 = 1. 

From Fig^ we notice that when s > 0, the optimized 
rate R over y increases with the degree of entanglement 
r and attains a maximum at some value r* > 0, so that 
the maximum is achieved by entangled input states and 
is greater than that at s — r = 0. 

Moreover, maximizing R with respect to both r and 
y allows to find the behavior of the transmission rate 
versus the number of channel uses n (see FigEJ) . Here 
we have considered s = 0,0.1,0.2, white, grey and black 
bars respectively. Unfortunately the evaluation of R with 
larger n requires a lot of numerical resources. Nonethe- 
less, from FigEI we can grasp that the transmission rate 



increases with the number of channel uses n, achieving 
an asymptotic value depending on s. 
R 

1.55 



1.50 




1 .40 



FIG. 2: The maximum of transmission rate R (over the input 
entanglement r and correlation coefficient y) as a function of 
the number of channel uses n for a channel with a degree of 
memory s = (white bars), s = 0.1 (grey bars), s — 0.2 
(black bars). The input mean number of photons is n = 2, 
while the number of added photons is N = |. 



This can also be understood by considering that the 
rate increases with noise correlations. Initially they 
rapidly grow till n reaches the relevant number of modes 
that effectively interact. In fact, each mode effectively 
interacts with a limited number of other modes which 
would be determined (besides the form of S) by the value 
of s (the smaller is s the smaller is the number of modes 
that effectively interact) . Thus the number of uses after 
which R becomes almost constant can be thought as the 
range of memory effects. 

In conclusion, we have provided a model to describe 
memory effects in a Gaussian channel with additive clas- 
sical noise for an arbitrary number of uses. The model 
reduces to that of Ref.Q for only two uses. From our 
numerical results we guess the existence of an asymp- 
totic value of the transmission rate after many uses of 
the channel that should lead to the full capacity of the 
channel. Hence, the 2-shot capacity would not suffice to 
characterize the channel performances, unless the mem- 
ory is very weak. Quite generally the minimum number 
of uses to consider would depend on the strength of the 
memory, i.e. the correlation length of the noise. 

Of course the results found are conditioned to to 
the validity of the conjecture that Gaussian inputs al- 
low to achieve the capacity of a Gaussian channel 
however, the presented model for memory effects is al- 
ways valid and could be extended to include attenua- 
tion/amplification besides additive noise. 
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